Can gravity distinguish between Dirac and Majorana neutrinos? 
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We show that spin-gravity interaction can distinguish between Dirac and Majorana neutrino 
wave packets propagating in a Lense-Thirring background. Using time-independent perturbation 
theory and gravitational phase to generate a perturbation Hamiltonian with spin-gravity coupling, 
we show that the associated matrix element for the Majorana neutrino differs significantly from its 
Dirac counterpart. This difference can be demonstrated through significant gravitational corrections 
to the neutrino oscillation length for a two-flavour system, as shown explicitly for SN1987A. 
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Introduction. -An unresolved question in the Standard 
Model is whether neutrinos exist as Dirac or Majorana 
particles [1, 2], where the latter are regarded as the more 
natural candidates [3] to exist in Nature. This is because 
the Majorana neutrino, as its own antiparticle, has only 
half the degrees of freedom compared to the Dirac neu- 
trino, and is deemed a more fundamental particle as a 
result. One known process to distinguish between the 
two types is neutrinoless double beta decay [1, 2], which 
only occurs for Majorana neutrinos. However, direct ob- 
servation of this phenomenon is at best inconclusive. As 
well, Dirac and Majorana neutrinos are potentially dis- 
tinguishable in magnetic fields, since they possess unique 
magnetic moment structures [1], though this feature be- 
comes relevant only when the fields are extremely strong. 
A recent discovery [4, 5] has shown that neutrinos un- 
dergo flavour oscillations while propagating in vacuum, 
inferring the existence of neutrino rest masses. However, 
the current theory of neutrino oscillations cannot distin- 
guish between a Dirac and Majorana neutrino because 
only its left-handed chiral projection is subject to this 
phenomenon [6], which is identical for both types. 

Although gravitational effects are usually neglected in 
particle physics, the fact that neutrinos are electrically 
neutral provides an opportunity to study their long-range 
behaviour in response to curved space-time. For example, 
it is possible that a massive neutrino's helicity state can 
be nipped due to explicit coupling between its spin and 
the background gravitational field. If such an interaction 
exists, then we have means to probe the intrinsic nature 
of neutrinos due to gravity, including the possibility to 
differentiate between Dirac and Majorana particles in a 
meaningful way. By treating the two neutrino types as 
massive wave packets propagating in a Lense-Thirring 
(LT) background [7], it is shown below that the Dirac 
and Majorana matrix elements differ significantly, along 
with their respective oscillation lengths. 

Dirac Hamiltonian in Curved Space- Time. -We be- 
gin with the four-dimensional covariant Dirac equation 
[i^(x)D^ - f ] ip(x) = of -2 signature (fj, = 0, 1, 2, 3) 



with neutrino mass m, with G = c = 1 units [8]. 
The curved space-time gamma matrices satisfying 
{'j^(x), J u (x)} — 2g t * l/ (x) are expressed in terms of 
Minkowski gamma matrices 7^ and orthonormal vier- 
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The LT metric [7] for x^ = (t, x, y, z) is 
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with M/r < 1 and MQR 2 /r 2 < 1, where r = 
^/x 2 + y 2 + z 2 , M and R are the mass and radius of the 
gravitational source, and is its rotational frequency. 
The corresponding Dirac Hamiltonian to leading order 
in M/r is 
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where a and (3 are the Dirac matrices, S J 
the a;- 7 -component of the spin angular momentum oper- 
ator, and L z is the orbital angular momentum operator 
in the z-direction. In spherical co-ordinates, the field 
point is r = (r,9,(p), defined in relation to a cartesian 
co-ordinate frame expressed by x^. The energy eigen- 
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where p = hko is the neutrino's momentum eigenvalue. 

Gravitational Phase. -By itself, (2) is insufficient to de- 
scribe a spin- 1/2 particle interaction with gravity. How- 
ever, for a weak field described by h^ v = g^ v — r\^ v <C 1, 
a gauge invariant gravitational phase 
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defined along some space-time trajectory z^ = 
(f, x', y', z') leads to a weak-field solution of the covariant 
Klein-Gordon equation [9-12], where p a and L a ^ are the 
generators of linear and orbital angular momentum for a 
free particle. Use of (4) according to the phase transfor- 
mation ip(x) — > exp (i^c/H) ip( x ) results in a new Dirac 
Hamiltonian H = Hq + H$ G , where 



H 9a = a • (V$ G ) + (V t $ G ) 



(5) 



is treated as a perturbation of Hq , reproducing all terms 
that have been observed, derived, or predicted to exist 
[13] for a spin- 1/2 particle in a gravitational field. 

Neutrino Wave Packets. -To formulate the wave packet 
description, we begin with a wavefunction composed of a 
linear superposition of plane waves [14] 
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where e lk ' r \U{k)) is a normalized solution of the free- 
particle Dirac equation, and 
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is a Gaussian function in momentum space, of width cr p 
and centroid fco> with fco = |feo|- The matrix element due 
to (5) is then 

(V>(r)|ff* G |V(r)) - ^3 J dV A 3 kA 3 k' 

x exp [i (k - k') ■ r'] (U(k')\H* G (r,r')\U(k)}, (8) 

where the integration is performed over all phase space 
in spherical co-ordinates, excluding the region occupied 
by the gravitational source. To evaluate (8) explicitly, 
we need to specify |£7(fc)). Assuming the Weyl represen- 
tation [2] for the gamma matrices, it is understood that 
the Dirac four-spinor is 
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where \vl(r)) is its left- (right)-handed chiral projection. 
In contrast, a Majorana four-spinor is identical to itself 
(up to some phase) under charge conjugation [2], where 
\ X C ) =±|X>- Then 
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After substituting (9) or (10) into (8), we proceed to 
evaluate the Dirac or Majorana matrix element with the 
Rayleigh plane wave expansion [15] in terms of spheri- 
cal Bcssel functions and spherical harmonics defined for 
both position and momentum space angles. It follows 
that the orthonormality conditions serve to truncate the 
series expansion, allowing a virtually exact evaluation of 
the matrix element. Generically, (8) is the sum of both 
spin-diagonal terms and spin-flip terms proportional to 
the Pauli spin matrices a 2 . Since the neutrino's spin 
quantization axis is parallel to the direction of propaga- 
tion, the helicity transition element is [14] 

(±|er|=p) = [cos 9 cos if ± i sin ip] x 

+ [cos 9 sin (f =F i cos ip] y — sin 9 z, (11) 

where |±) are the two-component spinors which define 
positive (negative) helicity for the neutrino. 

Dirac and Majorana Matrix Elements.-Aiter perform- 
ing a power series expansion of (8) with respect to 
?7i = m/(hk ) <C 1, we present the main formal results. 
For the Dirac neutrino, the gravitational phase-induced 
matrix element is 



(<A(r)|tf* G |?A(r)) Dirac = 

(Mo) [Co + C\ m + C 2 fh 2 
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where Cj and Dj are dimcnsionless functions of ho, R, r, 
and q = ko/c p , whose explicit expressions are shown in a 
much longer paper [16]. Some important details in (12) 
are as follows. First, Cj correspond to the spin-diagonal 
parts of (8) coupled to Mj r, while Dj refer to the spin-flip 
parts coupled to MQR 2 jr 2 . Second, the presence of sin 9 
clearly indicates that only terms with the ^-component of 
(11) survive the integration. The most obvious interpre- 
tation is that the gravitational source's rotation induces 
the helicity transition of the Dirac neutrino for propa- 
gation away from the axis of symmetry. That is, the 
off-diagonal metric terms in (1) resemble an inhomoge- 
neous magnetic field generating the spin-flip of a particle. 
Third, there are terms linear in fh present in (12) from 
the fact that the normalization coefficient in |J7(fc)) Dlrac 
is \J(E + m)/(2E), where E = yjp 2 + m 2 , with impor- 
tant consequences to follow. 
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(a) F 2 Dirac (SN1987A) 
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FIG. 1: F\ as a function of q due to the SN1987A gravitational 
source, for varying neutrino beam angle 9. Besides having 
opposite sign, J?j Dlrac is two orders of magnitude larger than 

pMaj. 
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FIG. 2: F2 as a function of q applied to SN1987A for varying 
9. The gravitational correction for F 2 Dlrac is larger than _F 2 Maj ' 
by three orders of magnitude. 



Evaluating (8) for the Majorana neutrino, we have 
(^ 1(2) (r)|^ G |^ (2) (r)} Ma 3- = 

(fifco)jy [Cq + C\ fa + C2 ?ri 2 ] 
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where the "1(2)" refers to the upper (lower) signs in (13), 
and Cjy and Djx are also dimensionlcss functions of fco, 
R, r, and q. The differences between (12) and (13) are 
quite striking. First, while both have the spin-diagonal 
terms Cj, the spin- flip terms have an overall factor of 
sin 8 sin ip. This corresponds to the y-component of the 
neutrino beam, a direct consequence of the Majorana 
neutrino's self-conjugation condition, as the charge con- 
jugation operation involves the presence of er y in its def- 
inition [2]. This result suggests a preferred direction or- 



thogonal to the source's axis of symmetry. However, since 
the LT metric is axisymmetric, this (^-dependence on (13) 
purely results from how we defined the co-ordinate sys- 
tem beforehand. This anisotropy can be removed by av- 
eraging over a complete cycle. Even so, if the source has 
a significant quadrupole moment to induce an azimuthal 
perturbation of the LT metric, then a resonance effect 
may be generated under suitable conditions. Second, 
the spin-flip parts of (13) are dependent on the x- and 
y-components of (11), as opposed to the z-component 
for the Dirac neutrino. Third, a spin-flip term still con- 
tributes to the Majorana matrix element in the limit as 
— > 0, while no such term survives for the Dirac coun- 
terpart. Though this seems counterintuitive, terms of 
this type are expected to be present because of the self- 
conjugate nature of Majorana neutrinos. 

Spin-Gravity Corrections to Neutrino Oscillation 
Length.-We now demonstrate how the matrix elements 
(12) and (13) lead to predicted gravitational corrections 
in the neutrino oscillation length, defined as L sc. = 



4 



2tt/ — E^^j [1, 2], where by convention we set 

> m\. To obtain , we use the Brillouin-Wigner 
(BW) method [17], instead of the more familiar Rayleigh- 
Schrodinger (RS) method [18]. From the BW approach 
applied to a second-order perturbation, we have 
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where the unperturbed energy Eq^ is described by (3), 
and (±\H& G |±) and (=p\H$ a |±) are the spin-diagonal 
and spin-flip components of the Dirac and Majorana ma- 
trix elements, as found in (12) and (13), respectively. 

The advantage of the BW method comes from knowing 
that it yields an exact expression for the total energy 
involving a two-level spin system, which is precisely what 
we have. Therefore, we can solve for E^ — E^ after 
averaging over the azimuthal angular dependence in the 
helicity transition term via 

|(T|^ G |±)| 2 - ^ fj \{t\H* g \±)\ 2 & V . (15) 
This leads to the expression 
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where F\ and F 2 are also dimcnsionlcss functions of fco, 
R, r, and q, and dependent on a complicated combination 
of terms coupled to M/r and MflR 2 /r 2 , whose analytic 
expressions are deferred to the forthcoming longer paper 
[16]. As a test case, we present F\ and F 2 as a function of 
q in Figures 1 and 2, respectively, for varying orientations 
of the neutrino beam angle 9, using data from SN1987A 
[19], where M w 1.4 M Q , R w 10 km, fl w 2.936 kHz, 
and r w 49 kpc. The neutrino wave packet is assumed to 
have a mean momentum of Hko = 1 MeV throughout. 

For all plots of Fj considered, the differences between 
positive and negative helicity are negligibly small. Also, 
the functions are insensitive to moderate or large changes 
in ho- Comparing Figures 1(a) and 1(b) for 9 — ir/2, we 
see that J\ Dirac and F\ Maj - have opposite sign, and \F^ iiac \ 
peaks near 3 x 10~ 2 , while the corresponding maximum 
i? 1 1 Maj - is roughly 6 x 10~ 4 . Regarding F 2 , Figures 2(a) 
and 2(b) for 9 = n/2 show that i^p irac has a maximum 
value near 0.5, while F^*' 1 ' peaks near 8 x 10~ 4 , about 
three orders of magnitude smaller. The fact that Fj are 
non-zero for the range of 10 -4 < q < 10 1 is consistent 
with reasonable choices of cr p for neutrinos produced in 
a neutron star, based on a mean free path calculation 
assuming known stellar data [14]. 

Conclusion. -This paper indicates that Dirac and Ma- 
jorana wave packets interact differently with gravity, with 



potential observational consequences. It is especially 
valuable to realize that the SN1987A values for F\, while 
small, are not negligible. This presents the interesting 
possibility that we can extract observational knowledge 
of the absolute mass difference m 2 — mi for a two-flavour 
oscillating system. Use of (16) can then lead to a deter- 
mination of the absolute neutrino masses m\ and m 2 by 
a parameter fit of q, m 2 — fh\, and m\ — m\ to precision 
measurements of the neutrino oscillation length, should 
this possibility become accessible in the future. 
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